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We study particle current and occupation distribution in a recently proposed model for coherent 
quantum transport. In this model a system connected to mesoscopic Fermi reservoirs (mesoreser- 
voir) is driven out of equilibrium by the action of superreservoirs with prescribed temperatures and 
chemical potentials described by a simple dissipative mechanism with the Lindblad equation. We 
compare exact (numerical) results for the non-equilibrium steady state particle current with theoret- 
ical expectations based on the Landauer formula and show that the model reproduce the behavior 
of coherent quantum systems in the expected parameter region. We also obtain the occupation 
distribution on the mesoreservoir in the non-equilibrium steady state and compare them with the 
occupation distribution on the leads in usual description of coherent quantum transport. 

PACS numbers: 05.60.Gg, 05.30.Fk, 03.65. Fd, 05.70.Ln 

I. INTRODUCTION 

The current of non-interacting particles flowing from a left to a right hand side reservoir trough a coherent mesoscopic 
conductor is usually described by the "Landauer picture" : Electrons in the left (right) reservoir that are Fermi 
distributed with chemical potential fi^ (rr) and inverse temperature Pl {/3r) can come close to the conductor and 
feed a scattering state that can transmit it to the right (left) reservoir. All possible dissipative processes such as 
thermalization occur in the reservoirs while the system formed by the conductor and the leads is assumed to be 
coherent. The probability of being transmitted is a property of the conductor connected to the leads, which is treated 
as a scattering system. In this picture, the probability that an outgoing electron comes back to the conductor before 
being thermalized is neglected, the contact is said to be reflectionless [IJ [2] . 

In the linear regime (small chemical potential difference) this scenario predict that G = ^%T{ep) is the conductance, 
where e is the electron charge, 2ttI% is Planks constant and T(ep) the total transmission probability evaluated at 
the Fermi level [T]. For simplicity we will consider spineless fermions but the generalization to the spinfull case is 
straightforward, in that case a factor two should be added to the conductance. For a perfect one-dimensional conductor 
T(cf) = 1 this formula predict a finite resistance. Instead Landauer [3] suggested that the conductance should be 
Go = i^T(ep) / R(tp) with R(ep) the reflection coefficient that vanish for a perfect conductor reestablishing the 
intuitive idea that an ideal conductor should have zero resistance. There was a long controversy regarding the validity 
of each formula that was finally settled with the realization of unavoidable contact resistance 1/G C , which is such that 
1/G C + 1/Go + 1/G C = 1/G. Associated to this contact resistance there is a potential drop [2111] taking place between 
the reservoir and the lead. Note that the resistance are added in series , therefore it is assumed that the contacts are 
not part of the coherent system. 

Recently several works have shown that the Landauer formula for the conductance or more generally the averaged 
particle current follows from a first principle description of the non-equilibrium steady state (NESS) of the non- 
interacting many-particle conductor-leads plus reservoirs system. In some particular limiting situations such as infinite 
reservoirs, this has been done rigorously [5] or using powerful tools such as non-equilibrium green function formalism 
that allows to treat the case where particles interact inside the system [6]. Most of these works focus their attention 
in the conductor and current but no attention is paid to the NESS properties of the leads or reservoirs, which are 
usually assumed to be in equilibrium. However, near the conductor deviations from equilibrium appear in reality as 
was recently observed in a cold-atom analog of a mesoscopic conductor [7]. 

In this paper we would like to discuss the relation between the current through a conductor and the non-equilibrium 
distributions on the leads in the framework of a recently proposed model [5] for a system plus mesoscopic reservoir 
that physically represent the leads that couple the conductor to the superreservoirs. Dissipation is included in our 
model and affects only the leads. The model is obtained by tracing out the superreservoir degrees of freedom in the 
unitary evolution equation for the density matrix of the full reservoir plus conductor system. At the price of several 
approximations such as the Born-Markov approximation (see, e.g., Ref. [51-U1)) a master equation in the Lindblad 
form for the reduced density operator of the conductor plus leads is obtained. The Lindblad dissipator encodes the 
effect of the Markovian superreservoir on the leads. 

The structure of the paper is as follows: In section |TI] we briefly present the model, in section |III| we analyze 
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the current flow and its agreement with Landauer formula and the occupation distribution on the mesoreservoirs in 



comparison with expectations based on the Landauer picture. Finally in section IV we present our conclusions 



II. DESCRIPTION OF THE MODEL 

As our "conductor" we consider both, a quantum dot and a one-dimensional quantum chain coupled at its left and 
right boundaries to finite reservoirs, which hereafter we call mesoreservoirs. The Hamiltonian of the total system can 
be written as H = Hs + Hl + Hr + V, where 

H s = UJc (1) 

is the Hamiltonian for the quantum dot or 

JV-l N 

h s = -J2 fah+i + ( h - c -)) + E U i c h ( 2 ) 
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is the Hamiltonian of the chain with tj the nearest neighbor hopping, Uj the onsite potential and Cj, cj the annihila- 
tion/creation operator for the spinless fermions on the site j of the chain (for the dot we do not need the subscript). 
The chain (or dot) interacts with the mesoreservoirs through the term 

K 

V = J2 ( v k a L c i + ^4i?c«) + (h.c), (3) 

k=l 

(in the dot case Ci and c„ should be replaced by c). The mesoreservoirs Hamiltonian are H a = X!fc=i £ fe a fe ct a fcQi here 
a = {L, R} denotes the left and right mesoreservoir. They corresponds to a finite number of spinless fermions with 
wave number &€{!,..., K} sharing the same spectrum with a constant density of states 9q in the band [E m i n , E n 
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described by = 9o(k — ko) and a.k tOL ,a\ are annihilation/creation operator of the a mesoreservoir. The system 
is coupled to the mesoreservoirs only at the extreme sites of the chain with coupling strength v? that we choose 
k- independent [T7] v% — v a . 

We assume that the density matrix of the conductor plus mesoreservoirs evolves according to the many-body 
Lindblad equation (hereafter we consider units where e = 1 and H = 1) 



= -i[H,p]+ ( 2L k, a , m pL{ am -{Ll am L k . a , m ,p}y (4) 
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where m G {I, 2} and Lk, a ,i = Vt(1 - fdfikjjaka, Lk,a,2 = \/lJjJk) a \ a are operators representing the coupling of 
the mesoreservoirs to the superreservoirs, / Q (e) = (e' 9cv ( e_AlQ ) + 1) _1 are Fermi distributions, with inverse temperatures 
P a and chemical potentials p a , and [•, •] and {•,•} denote the commutator and anti-commutator, respectively. The 
parameter 7 determines the strength of the coupling to the superreservoirs and to keep the model as simple as possible 
we take it constant. The form of the Lindblad dissipators is such that in the absence of coupling to the chain (i.e. 
v a = 0), when the mesoreservoir is only coupled to the superreservoir, the former is in an equilibrium state described 
by Fermi distribution [T^l [T3]. Physically the mesoreservoir mimic the leads (with a finite number of modes) that 
connect the macroscopic (infinite) reservoir to the system. This model has been justified starting from a Hamiltonian 
description of the baths in the case of the quantum dot |10j in the limit 7 — > 0, K — > 00 and we apply it also for short 
chains (N <S K). 

III. ANALYSIS AND PROPERTIES OF THE MODEL 

To analyze our model we use the formalism developed in |14) . There it is shown that the spectrum of the evolution 
superoperator is given in terms of the eigenvalues Sj (so-called rapidities) of a matrix X which in our case is given by 
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where Oi X j and Ej denote i x j zero matrix and j x j unit matrix, a y is the Pauli matrix, and H is a 
matrix which defines the quadratic form of the Hamiltonian, as H = wHd in terms of fermionic operators 
d T = {a 1L , ■ ■ ■ dKL>Cx, • • • , cjM,a 1R , ■ ■ ■ a KR }. 

The NESS average of a quadratic observable like d^di is given P3] in terms of the solution Z of the Lyapunov equation 
X T Z + ZX = Mi with M, = -|diag{rni L , • • • , m KL , 0i xN , mm, ■ ■ ■ ,m KR } (g> cr y and m ka = j{2F a (e k ) - 1} as 
follows: Consider the change of variables vi2j-i = dj + dp w%j = i(dj — G?j), the NESS average of the quadratic 
observable WjW k is determined by the matrix Z through the relation (wjWk) — Sj.k — 4&Zj,h- Wick's theorem can 
be used to obtain expectations of higher-order observables and in fact, the full probability distribution for these 
expectation values in some cases 15J. 

As shown in |14] if there is no rapidity with zero real part the NESS is unique. We verify this in our model and 
thus all the steady state properties that we describe bellow apply to the unique NESS of the system 



A. Particle current 



In this model, particle current is conserved for all parameters Heisemberg equation of motions for the occu- 
pations on the mesoreservoirs n a (e k ) = a\ a aka gives |16j 

^j_n a (e k ) = -j a (e k ) + V(n a (e k )) (5) 

with T>(n a (e k )) — —2~f(n a (e k ) — f a {e k )l) being interpreted as the current flowing from the a superreservoir to the 
k level of the a mesoreservoir and j a {<^k) — iv k a{a\ a c — c< a ka ) being the current flowing from the k level of the a 
mesoreservoir to the quantum dot or quantum chain considering c = c\ if a = L and c = c„ if a = R. 
In the dot j^n = J L + J R with J L = J2 k JL( e k) and J R = J2 k JR( e k)- And along the chain 

= J(_i - Ji l = l,...,n (6) 

where J ; = -it(c]c ;+1 - cj +1 c ; ) for I = 1, . . . , n - 1, J = J L and J n = -J R . 

In a steady state where all time derivatives vanishes we can see that ^ k {jh{^k)) + SfcO'ii( e fc)) = for the quatum 
dot and J2 k (jL{^k)) = (Jo) = ■■■ = (,/„) = —J2 k (jii( e k)) for the quantum chain. Thus the total average particle 
current is conserved and we can forget the chain sub-index and call the total average particle current (J). In our 
model generically (jhi^k)) + (ifl(efe)) 7^ 0. The equality is only found in a parameter region that we call the Landauer 
parameter region where transport is coherent as when collisions are elastic. When dissipation (7) is strong, the particle 
current is redistributed in different energy shells analogously as in inelastic scattering processes. Therefore, as we 
showed analytically in [TB], energy current is not generically conserved but total particle current is conserved. 



B. Particle current in the coherent regime: The Landauer Formula 

The Landauer formula [5] provides an almost explicit expression for the NESS (assumed to exist and being unique) 
average current in a coherent system as a function of the parameters of the system. In units where e = 1 and H = 1 
it reads 

(J) = il I de(f L (e)-f R (s))T(e), (7) 



2vr . 

where T(e) = tr[r l(e)G + (e)T R (e)G~ (e)] is the transmission probability written here in terms of 

G±{£) ^e-H s -4(e)-^ R (ey (8) 

the retarded and advanced Green function of the system connected to the leads and of —T a /2 the imaginary part of 
the self-energy S+. For the quantum dot with symmetric coupling at the left and right mesoreservoir [T5] 

Tie) = — ^ (9) 

[) ([/-e + 2A(e)) 2 +r(e) 2 y) 
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with Vqp(e) — r(e)/27T and 

i r md . = v> 

w 2vr J s-lj O 

In the so called wide-band limit E m i n <C U and C/ <C E max , A(u>) can be neglected. An explicit, but cumbersome 
expression of T(e) for the the quantum chain |15l ? ] can be obtained. We will use it for the numerical examples 
reported bellow but we do not reproduce it. 



(10) 



C. Mesoreservoir occupation distribution 

Let us consider Eq.([5]) for the left mesoreservoir in the NESS 

(j L (e k )) = -2 7 [(n L (s k )) - f L {s k )] (11) 

This equation nicely shows that current, a fingerprint of the non-equilibrium steady state, can flow only if the 
mesoreservoir distribution departs from the equilibrium distribution. We can rewrite this equation as 

(n L (s k )} = h{e k ) - ±{j L (s k )) (12) 
27 

which gives the NESS distribution on the left mesoreservoir if the current 0'l(£/c)) is known. When transport is 
coherent an analytic expression for {jhi^k)) can be given. In the general case, we resort to numerical computations. 
In fact, in the coherent case (ji(£fe)) = — (jR(£k)) and thus (jL(efe)) can be directly read from Landauer formula 
Eq.Q for the current 

PL(e k ) ]L (e k ) = ^-[h{e k ) - f R (e k )]. (13) 



From Eq.(12| and Eq.(13l we obtain 



and similarly 



n L (e k ) = f L {e k ) - _^g^ [/ £ ( £fe ) - f R (e k )} (14) 



n R (s k ) = f R (e k ) + J^^ fdeu) - f R (e k )] (15) 



which is one of the central results of this paper. For the sake of generality we consider in these formulas the generic 
case with different density of states in the left Phi^k) and right p R {e k ) mesoreservoirs but as we mentioned before in 
this work we consider them equal (p^ — Pr = l/$o)- We remark that Eqs.( 14|15 ) are valid if the current is given by 



Landauer formula, a situation that we discuss next with the help of numerical results. 



D. Numerical results 



The current: In Fig[T]we plot the total particle current (J) through the dot as a function of 7 for different values 
of K. For a fix value of K, we observe an initial linear growth of (J) with 7 consistent with the fact that for small 
7 the super- reservoirs are weakly coupled to the mesoreservoirs plus dot system and therefore the distribution 



differs from fa and according to Eq.(ll) the current grows linearly with 7. We also observe that if we consider a 
fix small value of 7 in this linear regime, (J) increases if we increase K (compare big blue dots with small red dots) 
but never goes beyond the value computed from Landauer formula Eq.|7| with T(e) from Eq. ([9]), indicated by the 
thin line in Fig[T] In fact we observe that when the current reaches the Landauer value it stays there in a plateau 
and then starts to decay with 7. Interestingly the curves obtained for different K superpose in this decaying regime. 
This provides interesting information about the strong couplin g re gime. For large 7 the mesoreservoir distribution 



approaches the corresponding Fermi distribution and from Eq.(ll) we can conclude that — fi ~ 0{~f~ / K). 
The plateau occurs at the maximum value of the current which is moreover given by the value predicted by Landauer 
formula. We observe that by choosing 7 and K such that 

2 1Pa = 2 7 /# » 1, (16) 
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FIG. 1: Particle current versus 7 for the quantum dot with U = 0, /3l = 5, /3r = 10 , y,h = 0.5, fiR — —0.5. The red dots 
are for K = 51, the blue dots are for K = 501 and the green dots are for 7 = #o/2 = (-E max — _E m i n )/(2if). The thin line 
indicates the value of the current pre dicted by Landauer formula for the same parameters (3l,Pr, Ml, P-r and vo- In this and 
al subsequent figures we take: vq = \J Oo / {2tv) and _E ma x = — -Emm = 5. 
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FIG. 2: Mesoreservoirs occupation distributions in the Landauer regime. The dots are computed numerically for the single 
dot model in the NESS a) Green dots: occupation distributions for the left mesoreservoir with K = 500. Continuous line is 
obtained from Eq.(14l with parameters chosen according to Eq.(16l. b) 
250. 



mesoreservoir with K 



Green dots: occupation distributions for the right 
Continuous line is obtained from Eq.(15l with parameters chosen according to Eq.(16l. 



the computed value of the current always fall on the Landauer value (green dots). Condition Eq.(16) defines the Lan- 
dauer parameter region mentioned above and is interpreted as equating the width (due to contact with environment) 
of each level of the mesoresevoir 27 with the level spacing 9 a . In other words, the density of states in the mesoreservoir 
is smooth. We have also check numerically (data not shown) that in this line (7 = #o/2) the current is conserved on 
shell i.e. jL^k) = —jn(£k), which is consistent with the fact that particle current is given by Landauer formula. 

Occupation: Now we analyze the mesoreservoir occupation distribution for the different regimes of (J) versus 7 we 
just discussed. 

• In the strong dissipation regime 7 > 6q/2, where the current start to decay as 7 increases (see Fig[T]), we expect 
that n a {tk) ~ /a(efe) as is verified in Fig. [3] 

• In the weak dissipation regime 7 < 9q/2, where the current grows linearly with 7 (see FigjTJ) , the super reservoir 
are almost decoupled from the mesoreservoir and the whole mesoreservoir-dot system is isolated. In that case, 
we expect from non-equilibrium Green function formalism that the occupation of the diagonal modes of the 
full Hamiltonian H = Hl + Hs + Hr + V is the average of left and right Fermi distribution. Thus for large 
mesoreservoirs (large K) we expect Ui(efe) = nji(ek) ~ (/i(efe) + //j(efe))/2 in the symmetric case considered, 
as verified in Fig[2J 



In the Landauer regime i.e. when 2j/9q 1 and the current is in the plateau, the non-equilibrium distribution 



is predicted by Eq.(14| with 2jpl = 1 i.e 



n L (e k ) = h(e k ) - ^[/^fe) - /flfab)] (17) 
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FIG. 3: Mesoreservoirs occupation distributions in the strong coupling limit 7 = 10. The dots are computed numerically for 
the single dot model in the NESS with K = 51. a) Red dots: occupation distributions for the left mesoreservoir. Continuous 
line under the dots is /z,(e) while the red line is the expected population in the Landauer regime obtained from Eq.(14l with 
parameters chosen according to Eq.(16l. b) Blue dots: occupation distributions for the right mesoreservoir. Continuous line 



under the dots is /ii(e) while the red line is the expected population in the Landauer regime obtained from Eq.(15l with 
parameters chosen according to Eq.(|16[). 
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FIG. 4: Green dots: mesoreservoirs occupation distribution in the weak coupling regime for the right mesoreservoir for 7 = 10~ J 
and K = 51. The continuous line under the dots is given by (/l + /r)/2, while the blue line is the expected population in the 
Landauer regime. The left distribution is identical to the right distribution and we do not show it. 



which is plotted in Fig|2]and perfectly agree with the numerical data. 

Note that the first two regimes are not physically relevant for our context here, which is to model the current and the 
leads in a typical non-equilibrium situation where the parameters 7 and K are not under control. 

We have done the same analysis for a chain. The same picture is observed for particle current as a function of 7 
and for the occupation of the mesoreservoirs, but since the transmission coefficient has an interesting structure as a 
function of energy, these structure is observed in the mesoreservoir distribution as illustrated in Fig(5] 

Physically the distributions n^^Sk) and nji(sk) on the mesoreservoirs correspond to the distribution on the leads, 
since mesoreservoirs model the leads that connect the superreservoirs to the system. The non-equilibrium distribution 
of the leads has been considered in [3] and in more detail in [2J. We compare both results bellow. 



E. Mesoreservoirs versus Fully coherent leads 



In the Landauer picture [3J 0] the contacts are reflect ionless thus the occupation distribution of electrons traveling 
to the right at the left lead is n\ = fa while at the right lead is = Tfi + RJr = Jr + T(fi — Jr). Analogously 
the distribution of left going electrons at the right lead is — Jr while at the left lead is = TJr + RJl = 
T(Jr — /l) + /l- From the difference of the distribution at the left and right we can evaluate the the potential drop 
across the scatterer (from left to right) for left and right moving electrons and obtain A/x c = (1 — T)(/X£ — ^r) at 
zero temperature. Accordingly we can define a conductance of the barrier Go which is related to the conductance G 
by Go = — jhr Su = Ti as originally suggested by Landauer. We can also compute the potential drop from the 
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FIG. 5: For the chain with N = 6, K = 101, tj — 1, Uj = 0, /3_l = Pr = 4, — — = 4a) current and b) mesoreservoirs 
occupation distribution in the Landauer regime. The red dots are numerical computations and the continuous line is from 
Eq.flTFf with T(e) computed in [15]. 



reservoir to the lead and the associated contact resistance 1/G C obtaining that 1/G C + 1 /Go + 1/G C = 1 /G. 

If we consider that in the leads for every right going state there is a left going states with the same energy then 
the occupation for the energy states on the left lead is fi, + (T/2)(/# — f£) which differ from in Eq.([T7"|). The 
difference is due to the following: In our model, the mesoreservoirs are subject to dissipation, so in the Landauer's 
view they are still part of the reservoir and therefore the potential drop associated to the difference between fr, and 
nr, is less than the potential drop associated to the contact resistance. In this sense mesoreservoirs are intermediate 
between reservoirs (described by Fermi distribution) and leads (that are not affected by dissipation) . Note that in the 
Landauer picture the system plus leads is considered to be coherent, thus contact resistance and system resistance 
are added incoherently. 



IV. CONCLUSIONS 



We have studied a model of coherent transport where the non-equilibrium distribution of mesoreservoirs can be 
analyzed. In this model Landauer prediction for the current holds under certain (and well understood) choices of the 
parameters and it is the maximum particle current of our model. We would like to point out that it may be possible 
to control dissipation in the experimental setting discussed in [7] , therefore other parameter regions (out of Landauer 
region) of our model could be relevant. 



The non-equilibrium distribution on the mesoreservoirs [Eqs. ( 14|15 ) in the parameter region defined by Eq.(16l] 



differs by numerical factors form the ones discussed in section II 



E| expected from Landauer's picture. As we com- 



mented above, in our model dissipation act on the mesoreservoirs and therefore they are not part of the coherent 
system. Indeed there is a potential drop between our distribution Eq. (14) and the one expected from Landauer's 
picture (see section III E I . In this sense with mesoreservoirs there is a smother potential drop between left reservoir 
and conductor. 
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